Astronomy & Astrophysics manuscript no. 1337bis 


© ESO 2009 


February 5, 2009 





Star formation around RCW 120, 
the perfect bubble * 

L. Deharveng^, A. Zavagno^ F. Schuller^, J. Caplan\ M. Pomares^ and C. De Breuck^ 

' Laboratoire d'Astrophysique de Marseille (UMR 6110 CNRS & Universite de Provence), 38 rue F. Joliot-Curie, 13388 Marseille 
Cedex 13, France, 

- Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany, 
^ European Southern Observatory, Karl-Schwarschild Strasse, 85748 Garching bei Miinchen, Germany 
e-mail: lise . deharvengOoamp . fr 



Received 12/11/2008; accepted 15/01/2009 

ABSTRACT 



o 
o 

(N 

Context. This study deals with the star formation triggered by H n regions. 

Aims. We wish to take advantage of the very simple morphology of RCW 120 - a perfect bubble - to understand the mechanisms 
^ triggering star formation around an H ii region and to establish what kind of stars are formed there. 

["T I Methods. We present 870 jjm observations of RCW 120, obtained with the APEX-LABOCA camera. These show the distribution 

of cold dust, and thus of neutral material. We use Spitzer-MlPS observations at 24 yum and 70 lum to detect the young stellar objects 
present in this region and to estimate their evolutionary stages. 

Results. A layer of dense neutral material surrounds the entire Hn region, having been swept up during the region's expansion. 
This layer has a mass greater than 2000 Mq and is fragmented, with massive fragments elongated along the ionization front (IF). 
We measured the 24 yum flux of 138 sources. Of these, 39 are Class 1 or flat-spectrum young stellar objects (YSOs) observed in the 
direction of the collected layer. We show that several triggering mechanisms are acting simultaneously in the swept-up shell, where 
• they form a second generation of stars. No massive YSOs are detected. However, a massive, compact 870 yum core lies adjacent to the 

IF. A 70 yum source with no 24 yum counterpart is detected at the same position. This source is a likely candidate for a Class YSO. 
Also at 24 yum, we detect a chain of about ten regularly spaced Class 1 or flat spectrum sources, parallel to the IF, in the direction of the 
JL most massive fragment. We suggest that the formation of these YSOs is the result of leans gravitational instabilities in the collected 

^ layer. Finally, the 870 yum emission, the 24 yum emission, and the Ha emission show the existence of an extended and partially ionized 

photodissociation region around RCW 120. This demonstrates the long-distance influence of the Hii region upon its surrounding 
medium. 
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1. Introduction IF. Many young stellar objects (YSOs) are present in the vicinity 

of the Hii region. However no massive object is observed, even 

g\ RCW 120 IS perhaps the most perfect of the many bubbles de- ^j^^^jj^j^ ^j^^ ^^^^ ^^^^-^^ fragment. ZA07's conclu- 

p tected m the Galactic plane by Spitzer at 8.0 ^m (Churchwell ^-^^ ^^^^ ^^jj^^^ ^ ^^jj ^^^^^^ ^^^^ j^^^ ^^^^^ 

^ EOOg bubble S7). This H ii region presents a simple morphol- ^^0, but that this region is too young for massive-star for- 

Q ogy, with a well defined ionization front (IF) separating the ion- ^^^-^^ ^^y.^ 1^^^ ^-^ gravitational instabihties along the col- 

ized and the surrounding neutral material. This optical H ii re- jgcted layer 

gion has the additional advantage of lying nearby, at a distance 

of 1 .34 kpc, and being isolated. Its distance is reliable as the pho- Since then, new observations have become available, allow- 

tometric and kinematic determinations are in good agreement to progress in our understanding of star formation tng- 

(Zavagno et al. [20021 hereafter ZA07). The J2000 coordinates S^red by RCW 120. These observations are: i) a deep map of the 

^Cl of its 08 centi-al exciting star are I = 348°2239, b = +0°4643, ^l^^t emission at 870 //m, obtained with the APEX-LABOCA 

a = 17'>12'"20i6, 6 = -38°29'26". This star is identified in camera during the science verification phase, giving improved 

ZAOT's fig 1 knowledge of the cold dust distribution and hence of the neu- 

Star formation associated with RCW 120 has been studied ^ral material associated with RCW 120; ii) Spitzer-MIPSGAL 

by ZA07, mainly using observations obtained at 1.2 mm, sup- fr^m^^ 24 /zm and 70 ^m, allowing a better detection of the 

plemented by 2MASS and Spitzer-GLIMPSE data. ZA07 have ^SOs possibly associated with RCW 120, and a better determi- 

shown the presence of a layer of collected neutral material sur- "^^i"" °^ *eir evolutionary stages. The 870 //m observation is 

rounding the ionized gas. This layer is fragmented, with massive presented in Sect. 2.1, and the 24 fim and 70 fim frames and 

fragments adjacent to the ionized gas, and elongated along the *eir photometi'ic reduction in Sect. 2.2. Sect 3 deals with the 

distribution of the neutral material around RCW 120. The nature 



Send ojfprint requests to: L. Deharveng of the YSOs associated with RCW 120 is discussed in Sect. 4. In 

* Table 1 is available in electronic form at the CDS via anony- Sect. 5 we try to understand the interactions between the central 

mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via ptp://cdsweb.u-| Hii region and the surrounding material, and how these interac- 

|strasbg.fr/cgi-bin/qcat?J/A-HA/| tions influence the star formation in this region. 
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Additional details about RCW 120, and various composite 
colour images, can be found in Deharveng & Zavagno (I2008I I. 

2. New observations 

2.1. APEX-LABOCA observations at 870 /jm 

The APEXQlABOCA camera was used at 870 yum (345 GHz) to 
observe the continuum emission of the cold dust associated with 
RCW 120. The Large Apex BOlometer CAmera, commissioned 
in May 2007, is a 295-pixel bolometer array developed by the 
Max-Planck-Institut fiir Radioastronomie (Siringo et al. I2007I I. 
The APEX beam size at 870 yum is 19'.'2. 

RCW 120 was observed in July 2007 as part of the science 
verification program. A 15' x 15' map was obtained in the spi- 
ral raster mode. The on-source observing time was less than two 
hours. During these observations the amount of precipitable wa- 
ter vapour was ~ 1 mm. 

The observations were reduced using the Bolometer array 
data Analysis package (BoA; Schuller et al. in preparation). The 
reduction steps involved in the processing of the data are: flux 
caUbration and opacity correction, deduced from skydip mea- 
surements and observations of primary and secondary calibra- 
tors; flagging of bad bolometers; correlated noise removal; de- 
spiking; low-frequency filtering. As a result of the correction 
for correlated noise, uniformly extended emission (on scales 
larger than 2!5), which mimics the variations of the sky emis- 
sion (skynoise), is filtered out. 

The calibration uncertainty in the final map is of the order of 
15%. The pointing uncertainty is ~ 4" (0.025 pc at the distance 
of RCW 120). We estimate the rms noise of the final image to be 
0.02 Jy/beam. 

2.2. Spitzer-MIPS observations at 24 fim and 70 fim 

Fig. |2.2| is a composite colour image of RCW 120, from the op- 
tical to the mid-IR. Blue is the Ha emission of the ionized gas, 
from SuperCOSMOS (Parker et al. 2005). Green is die 8.0 yum 
emission from Spitzer-GLIMPSE (Benjamin et al. l2003l l. At this 
wavelength we mainly see the emission of polycyclic aromatic 
hydrocarbons (PAHs). These big molecules are destroyed inside 
the H II region by the UV radiation of the exciting stars. They are 
excited in the photodissociation region (PDR) by the radiation 
leaking from the H ii region. Thus they are good tracers of ion- 
ization fronts. Red shows the emission at 24 fim, from Spitzer- 
MIPSGAL (Carey et al. 2005). We see both the extended emis- 
sion of the small dust grains present inside the ionized region and 
the PDR, and the emission of grains present in the envelopes of 
YSOs and evolved stars, hence the usefulness of the 24 yum im- 
ages for detecting YSOs. (The central 24 fim extended emission 
is saturated.) 

We have measured three Post BCD frames from the Spitzer- 
MIPSGAL survey at 24 fim. These cover a zone going from 
257?85 to 258°35 in right ascension and from -38^70 to -38^20 
in declination, encompassing the zone of influence of RCW 120. 
We used the DAOPHOT stellar photometry package with PSF 
fitting (Stetson 198J7). The main difficulty to overcome is the 
measurement of faint stars superimposed on a bright and very ir- 
regular background. For this we used an iterative process. A first 



' This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX), in run 078.F-9005(A). APEX is a col- 
laboration between the Max-Planck-Institut fiir Radioastronomie, the 
European Southern Observatory, and the Onsala Space Observatory. 



run of DAOPHOT, using the normal procedures, yields stellar 
positions and rough magnitudes. These stars are then subtracted 
from the original frame. The resulting image is smoothed with a 
median filter (window of 1 1x1 1 pixels, larger than the size of a 
star) to remove any vestiges of the stars, and obtain an image of 
the background emission alone. This background image is then 
subtracted from the original frame. On the resulting image the 
stars lie on a low brightness and more uniform background (the 
only remaining structures are those smaller than the filter win- 
dow). Further reduction is performed by the normal DAOPHOT 
procedures. This iterative process can be repeated several times 
until convergence is satisfactory. We check the result by eye; if 
the reduction is good no brightness peaks or holes should be seen 
at the positions of the stars after these stars have been subtracted 
from the original frame. 

The 24 fim frames overlap, providing us with double val- 
ues of the magnitudes of 27 sources. This allows us to es- 
timate the accuracy of the results: the rms magnitude differ- 
ence for these 27 pairs is 0.11 mag. The DAOPHOT PSF mag- 
nitudes are then converted into aperture magnitudes, follow- 
ing the "Quick Point Source Photometric Measurement" pro- 
cedure (http://ssc.spitzercaltech.edu/archanaly/quickphot.html), 
and using isolated bright stars. The resulting 24 fim magnitudes, 
obtained for 138 sources, are given in Table[T] 

We have also reduced three 70 fim Post BCD frames, from 
the Spitzer-MIPSGAL survey. Six sources have been measured 
manually, using aperture photometry. These measurements are 
also given in Table [T] 

We have used the interactive software sky atlas Aladirj^to 
superimpose the 2MASS catalogue on the GLIMPSE and 24 fim 
images. In the few cases where several GLIMPSE or 2MASS 
sources were observed in the direction of a 24 fim source, only 
the central one was considered. 

Table [T] which is available in electronic form at the CDS, 
presents our measured 24 fim magnitudes along with these ob- 
jects' 1.25 fim to 8.0 fim magnitudes taken from the 2MASS 
and Spitzer-GLIMPSE catalogues. Column 1 gives an identi- 
fication number Columns 2 and 3 are the coordinates accord- 
ing to the Spitzer-GLIMPSE catalogues. Columns 4 to 10 give 
the J, H, K, [3.6], [4.5], [5.8], and [8.0] magnitudes. These 
are extracted from the (highly rehable) GLIMPSE I Spring '07 
Catalog. When necessary they have been supplemented with val- 
ues from the (less reliable) GLIMPSE I Spring '07 Archive, 
and by the 2MASS All-Sky Point Source Catalog (a colon in- 
dicates an uncertain measurement). Columns 1 1 and 12 give our 
24 fim and 70 fim measurements. Five bright sources which have 
[24] < 1.38 mag (as indicated by a colon) are clearly saturated 
at 24 ftm. Column 13 gives our conclusions concerning the evo- 
lutionary stages of the sources (see Sect. 4). 



Table 1. Photometry of the 24 fim sources. Table 1 is available 
in electronic form at the CDS via anonymous ftp to cdsarc.u- 
strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi^ 
bin/qcat?J/A-HA/. 



^ http://aladin.u-strasbg.fr/ 
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Fig. 1. Colour composite image of RCW 120. Blue is Ha emission, green is the PAH emission at 8.0 jum and red is the 24 //m 
emission of small dust grains (see text). North is up and east is left. The field size is ~ 24' (N-S) x 25' (E-W). 



3. The distribution of neutral material around 
RCW 120 

The distribution of the 1 .2-mm dust emission was presented and 
discussed by ZA07 using SEST-SIMBA data (beam size 24", 
rms noise ~0.02 Jy/beam, exposure time 10 hr). Here we present, 
in Fig. |2] the distribution of the 870 jum dust emission, obtained 
with the APEX-LABOCA camera (beam size 19'.'2, rms noise 
~0.02 Jy/beam, exposure time < 2 hr). 

Emission at 870 /im is mainly the thermal continuum of cold 
dust. Contamination from the '^CO(3-2) line may occur, but this 
emission is expected to be faint (see the discussion in Schuller 
et al. in preparation). 

The 870 jum map clearly resembles the 1 .2-mm map: a layer 
of cold dust surrounds the ionized region; this layer is frag- 
mented, with condensations elongated along the ionization front, 
especially in the south. The higher resolution and higher signal- 
to-noise ratio at 870 jum reveal new structures that were only 
barely visible in the 1.2-mm map. This is the case for the ra- 
dial structures that are observed north-east of the ionized region, 
and for the filaments responsible for the absorption in the central 
parts of RCW 120 (the central absorbing filaments are conspic- 
uous on Fig. 2, right, blue channel). Also, small faint condensa- 



tions are now revealed, these often being associated with Class I 
sources detected thanks to the GLIMPSE and MIPSGAL sur- 
veys (Sect. 4). 

3.1. Mass determination 

Assuming that the 870 jum emission is uniquely thermal dust 
emission, we can derive the mass of the emitting material. 
According to Hildebrand jl983b the total (gas+dust) mass of a 
condensation is related to its flux density Syhy 

S 870 /jm 

M(gas+dust) = 100 — — , 

KSlOiim O870AimU dustj 

where D is the distance of the source, /fsvo/jm is the dust opacity 
per unit mass at 870 yum, and B87o^m(7^dust) is the Planck func- 
tion for a temperature Tdust- We have assumed a gas-to-dust ra- 
tio of 100. The dust temperature is taken as 20 K or 30 K, as 
often assumed for protostellar condensations (cf. Motte et al. 
2003, Johnstone et al. 12006b . Following Ossenkopf & Henning 
(1994) we adopt the value /(■870/jm = 1-8 cm^ g"'. The dust 
opacity depends strongly on the dust grains' coating (ice man- 
tles), compactness, coagulation evolution and the ambient den- 
sity (see their fig. 5a-c). This value agrees reasonably well with 
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Fig. 2. Left: 870 fim emission, observed with APEX-LABOCA. The contour corresponds to the 5 cr level of 0. 1 Jy/beam. The dust 
condensations are identified by a number used in ZA07 and in Table |2] Right: Colour composite image of RCW 120; the 870 jum 
emission appears in red, the Ha emission in turquoise. 



those adopted by various authors to derive masses from SCUBA 
850 yum emission maps (e.g. Kirk et al. 120061 Young et al. l2003l l. 

Table |2] lists the measured and derived properties obtained 
for several structures discussed in the text (Sect. 5). Column 1 
gives the structure's name. Condensations #1 to #8 are those 
identified in ZA07 (their fig. 4). Columns 2 and 3 give the emis- 
sion peak coordinates and column 4 gives the emission peak 
level. Columns 5 and 6 give the intensity value used to define 
the integration boundary and hence the flux density. The choice 
of the levels used to delineate the structures and measure their 
flux densities is somewhat arbitrary; they have been chosen in 
order to separate the different structures. Finally column 7 gives 
the derived masses depending on the adopted temperature, 20 or 
30 K; the lower mass corresponds to the higher dust tempera- 
ture. The masses are highly uncertain. As shown by column 7, 
the dust temperature uncertainty alone results in a mass uncer- 
tainty of a factor two. But the dust opacity also is uncertain by 
a factor of two (Henning, Michel, & Stognienko I2005I I. which 
implies an uncertainty of the same order for the mass. 

The masses we estimate from the 870 fim emission map are 
roughly twice those obtained by ZA07 from the 1.2-mm map. 
This discrepancy is partly due to the values adopted for the dust 
opacity, and mostly to the levels used to delineate the condensa- 
tions. The new values are more reliable, as the 870 fj.m map has 
a much better signal-to-noise ratio; also, the structures for which 
we estimate the masses are better defined in this paper. 

In order to derive the mass contained in the annular struc- 
ture that surrounds the ionized region (referred to as "Shell" in 
Table |2]l we integrated the flux over the surface enclosed by the 
5 cr level contour (corresponding to 0. 1 Jy/beam; see Fig.|2]i. The 
result is a total mass for the shell in the range 1200-2100 Mq. 

Condensation 1 contains a very bright core (emission peak 
of 13.3 Jy/beam) superimposed on a bright elongated plateau at 
a level ~ 2 Jy/beam. This core is resolved, with a deconvolved 



size of 14'.'6 X 7'.'7 (0.095 pc x0.05 pc; FWHM; see Sect. 5.3, 
condensation 1). 

In addition to the condensations listed in ZA07, a few other 
structures were measured, most of them associated with Class I 
sources. They are identified in Figs 9 to 14. 

3.2. Column density and extinction 

We have calculated the H2 column density, A^(H2), from the sur- 
face brightness F^-jo^^, using the formula 



K&lOfim B'ijOfim(Tdust) 2.3 Mn Qbeam ' 

where Fgvo^jm is expressed in Jy beam"', -Sgvo^m in Jy^ -^(1^2) is 
per square centimetre, the hydrogen atom mass mu is in grams 
and the beam solid angle ilbeam is in steradians. 

Adopting Qbeam - 9.817 10"^ sr corresponding to a beam of 
19'.'2 (FWHM), and assuming Tdust = 20 K, this gives 

A^(H2) = 3.136 10^2 F87o^„. 

From the classical relations A^(H -1- H2)/E{B - V) = 
5.8 10^' particles cm"^ mag"' (Bohhn et al. 1978) and Ay - 
3.1 EiB - V), we obtain Al/ = 5.34 10-^2 A^(H2). 

The bright core in condensation 1, with its emission peak of 
13.3 Jy/beam, is responsible for a visual extinction of 225 mag 
(assuming a dust temperature of 20 K). Similarly, all the ma- 
terial enclosed by the 5 cr contour (Fig. |2] 0.1 Jy/beam level) 
corresponds to a visual extinction greater than 1 .7 mag. Dust fil- 
aments - very conspicuous in front of the ionized region - are 
responsible for the absorption observed at optical wavelengths 
in the direction of the centre of RCW 120 (e.g. at Ha, in Fig.|2]i. 
The material observed at 870 fim in the direction of the ex- 
citing star of the Hii region has a small column density, cor- 
responding to a visual extinction < 1 mag. This is less than 
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Table 2. Mass estimates 
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the extinction, Ay = 4.65 mag, determined by Avedisova & 
Kondratenko ( |1984| l for the exciting star. The relatively large 
beam of APEX-LABOCA smoothes out small-scale local fea- 
tures of the dust emission. Also, large-scale structures are re- 
moved during the reduction (Sect. 2). The visual extinction is 
probably larger than estimated using APEX data, by some 3 - 
4 mag. 

The Spitzer emission maps at 8.0 fim and at 24 fim show, in 
absorption, several infrared dark clouds (IRDCs) in the vicinity 
of RCW 120. When we compare these maps with the dust emis- 
sion map at 870 yum. Fig. [3] we see that most of the cold dust 
condensations in the north have IRDCs as counterparts, whereas 
in the south the brighter cold dust condensations have no IRDC 
counterparts. This gives an indication of the respective positions 
of the dust condensations and the ionized gas: in the north, the 
dust condensations are in front of the PDR, whereas they are 
slightly behind the ionization front in the south. 

The visual extinction can be estimated from the 8.0 fim ab- 
sorption (cf. Vig et al. 120071 ), and compared with that found from 
the 870 jum emission. The Indebetouw et al. ( 2005 ) extinction 
law gives Ag.o^mMA: - 0.43; the Rieke & Lebofsky ( Il985b ex- 
tinction law gives Ax/Ay = 0.112. Thus Ay - 22.55 T8^„„ with 
Ts/im the optical depth at 8.0yum. If F and Fq are the 8.0 yum sur- 
face brightnesses, respectively in the direction of the absorbing 
cloud and in the direction of the surrounding background (as- 
sumed uniform), then Tg^,,, = ln(Fo/F). Of course, this method 
for deriving the extinction does not work if the absorbing struc- 
ture has too great an optical thickness; an Ay of 50 mag seems a 
reasonable maximum value. 

Fig. |4] allows us to compare these two determinations of 
the extinction in the direction of a zone situated north-east of 
RCW 120 (in the direction of condensation 5). For this region 
we have estimated Fq to be equal to 45 MJy/sr. The first con- 



tour, on each map, coiTesponds to a visual extinction of 1 .7 mag. 
These two maps are similar, with a few differences: 

- The region with Ay > 1.7 mag is more extended on the 
8.0 /im extinction map than on the 870 yum map. In the low 
extinction zones the extinction derived from the dust emis- 
sion at 870 fim is 1 - 2 mag less than that derived from the 
8.0 yum absorption. In the large extinction zones the two de- 
terminations of the extinction are in rather good agreement. 

- The 8.0 fim extinction map shows small-scale structures, 
mainly filaments and curls. In the 870 yum observations these 
structures have been smoothed out because of the larger 
beam. Note that many of the large extinction filaments are 
perpendicular to the ionization front; we will return to this 
point in Sect. 5.4. 

An IRDC is also observed in the 24 fim image, in the same 
direction. Because the surrounding background is far from uni- 
form, we did not try to estimate the extinction via the absorption 
at 24 jum. 

What can we learn about the extinction from the near-IR data 
(from the 2MASS survey)? This data is difficult to interpret in 
terms of extinction. Zones devoid of stars, hence directions of 
large extinction, are clearly present all around RCW 120 (ZA07, 
fig. 3), in the J, H and K maps. This is the case, for example, 
of the north-east large extinction zone shown in Fig. 3. Very 
few stars are detected towards these large-extinction filaments, 
whereas elsewhere numerous stars are observed, many of them 
displaying large extinction (up to 30 mag in Ay) in the J - H 
versus H - K diagram (not shown here). These are most prob- 
ably evolved stars situated in the distant background, possibly 
in the Galactic bulge (RCW 120 is less than 12° away from the 
Galactic centre). 
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Fig. 4. Visual extinction in the direction of the north-east condensation 5. The grey scale is the same for the two figures: white is for 
zero local extinction, black is for an extinction of 10 mag or larger Left: extinction estimated from the 870 yum emission. Right: 
extinction estimated from the 8.0 fim absorption. In each figure the contours correspond to an extinction of 1.7 mag and 5 mag. 




Fig. 3. Dust condensations observed at 870 fum and IR dark 
clouds in the vicinity of RCW 120. The red contours (0.1, 0.2, 
1.0, 2.5, and 5.0 Jy/beam) correspond to the 870 yum emission. 
The grey background is the Spitzer-GLIMPSE image at 8.0 fim, 
showing several IRDCs. The most massive condensations (in the 
south) are not IRDCs. 



4. Young stellar objects in the vicinity of RCW 120 

We wish to determine whether each source detected at 24 yum 
is a star - main-sequence or evolved - or is a YSO - Class I if 
its luminosity is dominated by an accreting envelope or Class II 
if its luminosity is dominated by a disk. (Robitaille et al. 120061 
prefer to describe YSOs as stage I or stage II sources.) 



In the absence of a spectroscopic signature, several indica- 
tors can be used to establish the evolutionary status of a YSO, 
depending on the wavelength range. 

- One indicator is the presence of a near-IR colour excess in 
the J - H versus H - K diagram. This excess of near-IR 
emission is generally attributed to hot dust associated with a 
disk (Lada & Adams [T992T l. 

- Another indicator often used with the Spitzer-GLIMPSE sur- 
vey is the position of the source in the Spitzer-IRAC colour- 
colour diagram, [3.6]-[4.5] versus [5.8]-[8.0]. In this dia- 
gram, unreddened stars (both main-sequence and evolved) 
have colours near (0,0), whereas Class I and Class II sources 
occupy different regions, as shown by Allen et al. (I2004I I and 
discussed by Whitney et al. C2003all2003bll2004b . 

- The slope of the spectral energy distribution (SED) between 
the near-IR and the mid-IR, for example from 2 //m to 10 
or 20 yum, is another indicator. This slope is defined as 
a - d(log(vFv))/d(log(v)) In the following we shall con- 
sider that Class I sources have an increasing SED, i.e. a > 0, 
Class II sources have a decreasing SED with -1.6 < a < 0, 
and Class III sources have a < -1.6. One can use for ex- 
ample the K versus /r-[24] diagram (Rebull et al. 200'Ql. 
Considering sources corrected for the interstellar extinction. 
Class I objects have K — [24] > 7.54, and Class II objects 
have K - [24] between 3.37 and 7.54. Class III and stars, 
which we do not attempt to separate, have K - [24] < 3.37. 
(Ordinary stellar photospheres have K - [24] ~ 0.) 

We shall compare the information given by these three indi- 
cators to discuss the evolutionary stages of the sources observed 
at 24 jum in the vicinity of RCW 120, and listed in Table [T| We 
shall show the limitations of these indicators and try new ones. 
We have been greatly helped in this study by the comprehensive 
discussion of Robitaille et al. (i2006i) . 

Fig.[5]presents the K versus A"- [24] diagram obtained for the 
118 sources that have K and [24] magnitude measurements. We 



^ Rebull et al. ( 120071 also identify "flat spectrum" sources, which are 
intermediate between Class I and Class II, and correspond to -0.3 < 
a < -1-0.3 (6.75 < A' - [24] < 8.31). The domain of flat spectrum 
sources is indicated in Fig. 4 
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Fig. 5. The K versus K - [24] diagram. Dilferent classes of 
sources are separated by dashed lines. The dotted lines corre- 
spond to the detection limits of 14.5 mag in K and 5.5 mag in 
[24]. (Note that these detection limits are not completeness lim- 
its.) The black arrow is a reddening line corresponding to a vi- 
sual extinction of 30 mag. Class I sources appear as large red 
dots. Class II sources as medium green dots, and all the other 
sources as small black dots. The same symbols will be used for 
the same sources in Figs|6]|7] and[8](left). 



are tempted to identify 24 Class I sources (red symbols), and 5 1 
Class II sources (green symbols) in this diagram. However, note 
the strong effect of interstellar extinction on the ^^-[24] colour, 
and the suspicious concentration of sources near the separation 
line between Class II and Class III objects. We shall return to this 
point later in this section. (Following Robitaille et al. 120061 we 
here refer - unconventionally - to the extinction as "interstellar" 
when it is due to dust located outside the envelope and/or disk of 
a source.) 

Fig. [6] presents the [3.6]-[4.5] versus [5.8]-[8.0] diagram 
for the 103 sources that have been measured in the four IRAC 
bands. We use the same symbols as in Fig. |5] for the Class I 
and Class II sources already identified with the first indicator 
(K versus K - [24]). Fig. |6] shows an almost perfect agreement 
between the two indicators for Class I sources. Eleven new very 
red sources (black squares), which have no measurable K magni- 
tude and hence are absent in Fig.|5] are identified as Class I. The 
situation is less satisfactory for the candidate Class II sources 
(according to our first indicator). A number of these sources lie 
inside the ellipse of stars. We shall return to this point later in 
this section. 

Fig.|7]is the J - H versus H - K diagram for the 1 12 sources 
with measured 7, H, and K magnitudes. All the Class I sources 
(identified by our first indicator) exhibit near-IR excesses, and 
they are nearly the only ones to do so. (The relative amounts 
of the near-IR excesses is uncertain, as almost all these sources 
have uncertain 2MASS magnitudes.) The fact that most of the 
sources previously classified as Class II do not have near-IR ex- 
cesses is surprising, since in models of Lada & Adams (1992), 
the hot dust situated in the disk emits in the near-IR. One possi- 
ble explanation is that the slope of the extinction curve is steeper 
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Fig. 6. The [3.6] - [4.5] versus [5.8] - [8.0] diagram. The black 
arrow is a reddening line coiTesponding to a visual extinction of 
30 mag. Each source is represented by the same symbol as in 
Fig.|5] The black squares are for sources lacking K magnitudes. 
The locations, according to Allen et al. (2004), of sources dom- 
inated by their photospheres (stars and Class III objects) and of 
Class II sources are indicated by an ellipse and a rectangle, re- 
spectively. 



than that adopted in Fig. [7] For this figure we have used the in- 
terstellar extinction law of Indebetouw et al. (12005 l l. According 
to these authors. Ay /A = 2.50 + 0.15 wAAhIAk = 1.55 ±0.08, 
giving a slope (Ay - Ah) I (Ah - Ak) - 1.73 ± 0.67. This range 
of values is illustrated in Fig.|7] showing that a small increase of 
the slope would result in a near-IR excess for numerous Class II 
sources. 

The sources with /r-[24] < 3.37 are Class III objects or stars 
- in the later case mainly evolved stars. Main-sequence stars at 
the distance of RCW 120 or in the background (thus possibly 
reddened) should have K > 5.67 mag (Martins & Plez 2006). 
Thus all the sources brighter than 5.67 in K and with K - [24] < 
3.37 are possibly giants or supergiants. This is compatible with 
their positions in the J - H versus H - K diagram. One of these 
sources, #83, is very bright in K and lies in the direction of the 
ionized gas. This star is HD 155275, an MOIII star (Houk [T982l 
it is identified in Deharveng & Zavagno'2008). 

Many sources may be affected by a large interstellar extinc- 
tion. Fig. |7] for example, shows that many stars (black dots) or 
candidate Class II sources (green dots) may be affected by an 
(external) extinction between 20 and 40 mag. This extinction 
may severely affect the classification of the sources. We suggest 
that some stars or Class III objects, affected by strong extinction, 
are mistaken for Class II objects by our first indicator. Our first 
indicator shows a large concentration of stars and Class III and 
Class II objects at K- [24] ~ 3.37. All the Class II sources which 
have 3.37 < K - [24] < 4.0 lie inside the ellipse (showing the 
location of stars) in Fig. [6] F urthermore, according to their 7-// 
and H — K colours (Fig. |7]i, all these sources may be evolved 
stars. Thus most of these stars are probably evolved background 
stars and not Class II sources. 

And the situation is perhaps not better for Class I sources. 
Fig. [8] presents the [24] versus [8]-[24] diagram, very simi- 
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Fig. 7. The J-H versus H-K diagram. The main sequence (MS) 
is drawn with a solid curve; the colours are from Martins & Plez 
(2006 ) for O stars, and from Tokunaga (120001) for later spectral 
types. The reddening lines (dashed lines), for 05V and M2V 
stars, are for a visual extinction of 50 mag (extinction law of 
Indebetouw et al.'2005 ). The symbols are the same as in Fig.|5] 
All the sources situated below the lower (O-star) reddening line 
exhibit a near-IR excess. The sources situated above the upper 
reddening line are possibly evolved stars. The solid-line arrows 
are the reddening lines corresponding to a visual extinction of 
30 mag and compatible with the extinction law of Indebetouw et 
al. ([2005. see text). 



lar to the K versus /r-[24] diagram of Fig. [5j with the disad- 
vantage of a smaller wavelength range but the advantage of a 
weaker influence of the extinction. In this diagram unreddened 
Class I sources have [8] - [24] > 3.58, and Class II sources 
have 1.67 < [8] - [24] < 3.58. We have kept the same sym- 
bols as given by our first indicator. Fig.[8](left) shows that some 
sources previously classified as Class I now appear as Class II, 
and that some sources previously classified as Class II are prob- 
ably Class III sources or stars. 

We shall now try to estimate the possible influence of the 
interstellar extinction on our classifications. We use the 870 fim 
emission to derive a maximum value of this external extinction; 
we assume that the sources lie in the centres of the dust conden- 
sations observed at 870 //m in their direction. The extinction is 
then given by 



Av(external) 



16.75 Fa- 



70 /i m 



+ 3, 



where F870/jm is expressed in Jy/beam and the 3 is to correct for 
a 3 mag underestimate of the extinction using the 870 fj.m map 
(Sect. 3.2). Fig. [8] (right) gives an idea of the maximum exter- 
nal extinction possibly affecting the sources. It shows that only a 
few sources are possibly misclassified for this reason. These are 
mainly sources observed in the direction of the bright conden- 
sation 1 (#41, #44, #46, #47, #49, #54) and in condensation 6 
(#106, #110) (Sect. 5.3). These sources may be flat spectrum 
sources. For the other sources, the discrepancy between the clas- 
sifications given by different indicators is most probably due to 



difficult photometric measurements. It mainly affects sources ly- 
ing near the separation limits of the different classes. 

Table [T] gives our conclusions concerning the evolutionary 
stages of the sources. As Class I we include 21 sources satis- 
fying all the indicators (K versus K - [24], [3.6] - [4.5] versus 
[5.8] - [8.0], [24] versus [8]-[24]) and not affected by too much 
external extinction. We consider as intermediate Class I-Class II 
18 sources which are Class I or Class II depending on the in- 
dicator (identified by "I-II?" in column 13 of Table 1), or are 
possibly flat spectrum sources mistaken for Class I because of a 
large external extinction ("I-flat?" in Table 1)). We consider as 
Class II 29 sources satisfying all the indicators, and as interme- 
diate Class Il-Class III 9 sources which are Class II or Class III 
depending on the indicators ("II-III?" in Tablel). Note that in 
Table [T] we have classified as "stars" all the sources dominated 
by their photosphere, and which we are unable to distinguish: 
Class III YSOs, MS or evolved stars. 

The spatial distribution of the YSOs is presented on Fig. [9] 
and is discussed in Sect. 5.3. 



5. Discussion 

5.1. Dust condensations and infrared dark clouds 

RCW 120 demonstrates in an exemplary way that massive dust 
condensations are not all infrared dark clouds. A dense neu- 
tral condensation needs to be seen towards a bright background 
emission to appear as an IRDC. But another necessary condition 
is to not be located behind bright foreground emission. In the 
case of RCW 120 the structures that are the most massive (con- 
densations 1 and 2) and the densest (the core in condensation 1, 
responsible for an extinction ~ 225 mag) are not IRDCs: they 
are masked by the bright emission of the adjacent foreground 
PDR. Only their external parts (those far from the IF) appear as 
faint absorption features (Fig. [3]). 



5.2. Young stellar objects and AGB stars 



Robitaille et al. ( 120081 ) present a catalogue of intrinsically red 
sources, observed by Spitzer, in the Galactic mid-plane. These 
sources consist mostly of high- and intermediate-mass YSOs 
and asymptotic giant branch (AGB) stars. As shown by these au- 
thors, it is very difficult to distinguish between these two types. 

In the following we compare the sources detected at 24 //m 
in the vicinity of RCW 120 with the red sources selected 
by Robitaille et al. If we apply their selection criteria, i.e. 
[4.5]-[8.0]> 1, 13.89 > [4.5] > 6.5, and 9.52 > [8.0] > 4.01, to 
the 112 sources with [4.5] and [8.0] magnitudes in Table 1, we 
find 39 sources; 26 of them are Class I or flat spectrum sources, 
1 1 are Class II, and two more are uncertain Class Il-Class III. 

Robitaille et al. estimate that ~19 percent of the red sources 
may be AGBs. If true, this means that up to seven of our 39 red 
sources may be AGBs and not YSOs! It is very difficult to distin- 
guish between these two categories if we only have photometric 
measurements up to 24 //m. Robitaille et al. show that YSOs are 
generally redder than AGBs between 8 //m and 24 fim, and pro- 
pose to use the criterion [8.0]-[24]> 2.5 to identify YSOs. Using 
this criterion we find that 32 of our 39 red sources are proba- 
bly YSOs (the seven sources with [8.0]-[24]< 2.5 are Class II or 
uncertain Class II). Thus we must keep in mind in the follow- 
ing that a few of our sources identified as YSOs, and especially 
Class II YSOs, are possibly AGB stars and are therefore evolved. 
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Fig. 8. The [24] versus [8]-[24] diagram, with the different classes separated by dashed lines. Class I sources have [8]-[24]> 3.58, 
Class II sources have 1.67 < [8]-[24] < 3.58. The reddening vectors are for a visual extinction of 30 mag. Left: The red and green 
dots correspond to Class I and Class II sources, according to our first indicator. Right: The maximum external extinction possibly 
affecting the sources is indicated as follows: red squares are for sources with Ay(external) > 30 mag, green squares for sources 
with 20 mag < Av(external) < 30 mag, blue squares for sources with 5 mag < Av(extemal) < 10 mag (no sources with 10 mag 
< Av(external) < 20 mag are present in this diagram). 



5.3. Large scale distribution of tlie YSOs 

We have measured 138 sources around RCW 120 emitting at 
24 fim. Three are unclassified. Twenty-one are Class I YSOs, and 
eighteen are possible Class I or flat spectrum sources. Twenty- 
nine sources are Class II YSOs, and nine sources are more uncer- 
tain Class II sources. Thus 49 to 56 percent of the 24 jim sources 
are young stellar objects, 28 percent of the 24 fim sources being 
Class I or flat spectrum sources. 

Fig.[9]shows the large-scale distribution of the YSOs possi- 
bly associated with RCW 120. It shows a large number of Class I 
or flat spectrum sources (red symbols), all but five associated 
with the shell of dense neutral material surrounding RCW 120. 
Thus even if we have no direct proof that any given YSO is as- 
sociated with RCW 120, this general picture demonstrates that 
most of them are associated with the H ii region, and that their 
formation has been triggered, one way or another, by the H ii re- 
gion. On the other hand, the Class II sources are more widely 
scattered over the field. Their association with RCW 120 is less 
evident, and this raises questions about their origin. 

We note the large concentration of Class I or flat spectrum 
YSOs in the direction of condensation 1, the most massive. 

5.4. Stellar content of the condensations 

In the following we discuss the YSO content of a few condensa- 
tions. 

• Condensation 1: This is the most massive of the condensa- 
tions, in the range 460-800 Mq. It contains a very massive core. 
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of 140-250 Mq, and several other bright substructures. Fig 
shows this condensation and the numerous 24 jim sources oF 
served in its direction or vicinity. 

The bright 870 jum core has a size (beam-corrected FWHM) 
of 14'.'6 X 7'.'7 (0.095 pc x 0.05 pc), along the IF and per- 



pendicular to it, respectively - thus slightly elongated along 
the IF. Its size and peak flux indicate a mean density of a few 
times 10^ molecules cm"-'. No bright 24 fj.m source is observed 
in the direction of the 870 fim emission peak, but three faint 
sources (including #47) may be associated with it. The Spitzer- 



MIPSGAL image at 70 yum in Fig. 1 1 shows an enhancement of 



the emission in the direction of the core, suggesting that a very 
deeply embedded object is present here. We measured a flux of 
130 Jy for this object, but this is very uncertain due to the un- 
derlying emission of the PDR and the low quality of the 70 //m 
image. 

According to our indicators, many 24 fim sources appear to 
be Class I YSOs. Among these are #38, which is observed in the 
direction of a spherical dust clump (b in Table |2] and Fig. [TOji; 
#58, also observed on the border of clump a2 in condensation 7; 
and #61, observed in the direction of the ionized gas (but most 
probably situated behind the H ii region). Sources #36 and #37, 
also located in projection close to a condensation, are possibly 
in an intermediate stage between Class I and Class II. The stage 
of #39 depends on the indicators used, but is probably a Class II. 
Source #32 has only [4.5] and [24] magnitudes; the slope of its 
SED suggests a Class I source. 

Condensation 1 contains a most remarkable feature: a chain 
of eleven 24 yum sources. They are aligned parallel to the main 
IF and are regularly spaced, separated by some 0.1 pc. They are 
deeply embedded. Sources #42 and #53 are Class I. Sources 
#41, #44, #46, #47, #49, and #54 lie neai- the separation be- 
tween Class I and Class II and may be flat spectrum sources as 
they may be affected by a large external extinction (cf. Fig. |8]l. 
Measurements lack for #40, #48, and #52; the slopes of their 
SEDs suggest Class I sources. 

We have fitted the SEDs of a few sources using 
the Web-based SED fitting-tool of Robitaille et al. (l2007l 
http://caravan.astro.wisc.edU/protostars/J . The parameters of the 
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Fig. 9. Distribution of the Class I (red circles), intermediate 
Class 1-Class 11 or flat spectrum (red crosses). Class 11 (green 
circles), and possible Class 11-Class 111 sources (green crosses). 
The grey-scale background is an unsharp-masked 24 jum image, 
where the extended dust emission has been subtracted to en- 
hance the stellar component. The contours represent the 870 fim 
emission level of 0.1 Jy/beam (5 cr). 



best model, as well as a range of possible values for some param- 
eters, are given in Table |3] YSO #38 is a stage 1 source, with a 
central object of a few solar masses; condensation b is probably 
not the source's cold envelope since its flux density at 870 fim is 
too large. Source #61 is also in stage 1, with a central object of 
a few tenths of a solar mass. For sources #53 and #49, most of 
the models indicate that they are stage 1 sources around central 
objects of about one solar mass, and with massive disks; but a 
few of the good models (those with^^ "A'best - show central 
objects with a higher mass, no accreting envelopes, and massive 
disks, hence corresponding to stage 11 sources. 

The chain of sources aligned along the IF, and regularly 
spaced, indicates that star formation has occurred in the col- 
lected layer, probably resulting from gravitational instabilities 
along a dense filament. At a temperature of 20 K and a density 
of 10^ molecules per cm^, gravitational instabilities form cores 
separated by the Jeans length Aj ~ 0.1 pc, and with the Jeans 
mass Mj ~ 1.5 Mq. This is probably what has occurred in- 
side condensation 1 . On the other hand the formation of the very 
massive and dense core may result from the collect and collapse 
process, i.e. from large-scale gravitational instabilities along the 
surface of the collected shell. Here again triggered star forma- 
tion is probably at work, as shown by the presence of a very 
deeply embedded source - hence a Class YSO - only detected 
at 70 fj.m. 

• Condensation 2: This lies at the south-east border of 
RCW 120, adjacent to the main IF. Its mass is in the range 110- 
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Fig. 10. Condensations 1 and 7: the background image is an un- 
sharp masked image obtained from the Spitzer-MIPSGAL frame 
at 24 yum. This image, free of the diffuse extended emission, en- 
hances all the stellar sources. The contours coiTespond to the 
870 jum emission. The contour levels are 0.2 Jy/beam (used to 
estimate the masses of condensations 1 and 7), 0.4, 0.8, 1.75, 
2.5 (used to estimate the mass of the core in condensation 1), 5, 
and 8 Jy/beam. 



Fig. 11. Colour composite image of the southern IF in the 
mid-lR. The Spitzer-MlPS emission at 70 yum appears in red. 
Turquoise shows an unshaip-masked image at 24 jum, obtained 
from Spitzer-MIPSGAL. The red lines point to the 70 //m source 
detected in the exact direction of the dense core in condensa- 
tion 1. 



290 Mq (not including the elongated structure which we call 
condensation 2bis in Fig. 12 1. Its peak column density coiTe- 

' 32 mag. The YSO #50 lies at less 



sponds to a visual extinction 
than 0.03 pc (in projection) from the peak direction. 

This source has no near-lR emission (in the 2MASS sur- 
vey), and no measured 8.0 fj.m flux, hence does not appear in 
our diagrams. However the slope of its SED between 3.6 jim 
and 24 jim is between 2.0 (without extinction correction) and 1 .3 
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(correcting for a visual extinction of 32 mag). Thus source #50 
is a Class I object. It lies very close to the IF, which is probably 
the reason why we do not detect it at 70 //m. 

Source #51 lies at the top of a "finger", a structure of the 
IF protruding inside the ionized gas, and pointing towards the 
exciting star This finger, which is seen south-west of condensa- 
tion 2, is an elongatd dust filament. Another such structure lies 
nearby, which we have called condensation 2bis (Fig. 12 1. ZA07 



proposed that these structures result from dynamical instabilities 
of the expanding IF, such as these simulated by Garcia-Segura 
& Franco ( 1996 ). These structures are found towards a gap be- 
tween condensations 1 and 2, through which the ionized gas tries 
to flow away from the central H n region (a region thus dynami- 
cally unstable). An alternative origin is possible. These two "fin- 
gers" are reminescent of the "pillars" in M16, with Class I & II 
YSOs located at their tops (Sugitani et al. 120021 Indebetouw et 
al. l2007l and references therein). The pillars in M16 clearly result 
from the interaction of NGC 66irs OB stars' UV radiation with 
pre-existing dense clumps. However, as shown by the molecular 
observations of White et al. ( 119991 ), most of the mass of the pil- 
lars is located in their heads rather than their tails, whereas the 
contrary is observed in RCW 120's fingers (at least in conden- 
sation 2bis). Thus the origin of the fingers and of YSO #51 is 
uncertain; higher angular resolution observations are needed. 

Our indicators show that YSO #5 1 is possibly a Class I ob- 
ject. Its SED, modelled by the Robitaille et al. fitting tool, con- 
firms this conclusion. The best model indicates a central object 
of 1.3 Mq with a massive disk but, however, with an accret- 
ing envelope, hence a stage I-stage II object. Several parameters 
are rather uncertain, but the disk is massive in all the models 
(Table|3]l. 

Finally the other YSO of the field, #59, is clearly a rather 
massive stage Il-stage III source. It is well fitted by a central 
object of 6.6 Mq, a small disk, and no envelope (Table [3]l. Its 
association with RCW 120 is uncertain. 

• Condensation 4: This is elongated along the IF but pro- 
trudes inside the ionized region. Its mass is between 50 and 
90 Mq. Fig. 13 shows that several 24 fim sources are detected 



in its vicinity, but that none is observed in the exact direction 
of the 870 yum emission peak. However all three sources (#64, 
#67, and #69) lie on the border of the condensation facing the 
ionized gas and have near-IR counterparts. Our indicators show 
that #64 is intermediate between Class I and Class II, and that 
#67 is a Class II. Source #69 has no measurement at 8.0 jim, so 
it is not in our diagrams, but the slope of its SED between 3.6 fim 
and 24 yum lies between 0.4 (uncorrected for extinction) and 0. 1 
(assuming Ay = 14 mag, corresponding to the 870 fim column 
density in its direction); thus source #69 is probably intermediate 
between Class I and Class II. 

We have used the Robitaille et al. SED fitting tool to model 
source #67. The result confirms our previous conclusion. The 
parameters of the best model are given in Table [3] This source 
contains a central object of mass ~ 5 Mq. Many of the good 
models have no accretion from an envelope. Thus this source is 
most probably a stage II, but the mass of the disk is very uncer- 
tain. 

Two regions of extended emission at 8.0 //m also lie on the 
border of condensation 4. ZA07 suggested that they were PDRs 
around Herbig Ae/Be stars. Emission from these PDRs is also 
observed at 24 fj.m. Their central stars have been observed by 
2MASS; their near-IR photometry suggests a B4V star affected 
by a visual extinction of 7.8 mag (object A in Fig. 13 or object 1 
in ZA07), and a B7V star with an extinction of 12.2 mag (object 
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Fig. 12. Condensation 2: the background image is the unsharp 
masked image at 24 yum; the contours correspond to the 870 //m 
emission. The contour levels go from 0.2 Jy/beam (used to es- 
timate the masses of the condensations) to 1.8 Jy/beam, by step 
of 0.2 Jy/beam. The dashed lines show the north-east and south- 
west limits adopted to delimit condensation 2. The nearby con- 
densation 2bis is identified. 



B in Fig. 13 or object 2 in ZA07). These stars are probably the 
most massive of the sources associated with condensation 4. 

ZA07 proposed that condensation 4 was a pre-existing dense 
clump which became surrounded by ionized gas during the ex- 
pansion of the H II region. Star formation is observed at its bor- 
ders and not at its centre. Thus star formation is probably not 
due to the clump implosion (which is predicted by the radiation- 
driven implosion models; Lefloch & Lazareff 119941 Kessel- 
Deynet & Burkert '2003] ). 

Fig. 13 shows the presence of a small spherical clump north 
of condensation 4, with a mass in the range 5-9 Mq. A bright 
24 yum YSO, #76, coincides with this clump to within the uncer- 
tainties (about 0.015 pc in projection). The evolutionary stage of 
#76 is uncertain: our indicators point to Class I or Class II. But 
most of the SED models have no accreting envelope, suggesting 
stage II. 

• Condensation 5: This is a group of structures emitting at 
870 fim on the north-east border of RCW 120 (Fig.[l4|. Its mass 
is in the range 130-225 Mq. It contains several small spherical 
dust clumps in the direction of which YSOs are observed: clump 
a, of mass 15-26 Mq, possibly associated with #107; clumps bl 
+ b2, of total mass 11-19 Mq, possibly associated with #105 
and #103. Clump c, of mass 6-1 1 Mq, possibly associated with 
#88, lies outside condensation 5. 

Numerous 24 fim sources lie in the field of Fig. [14] 
According to our indicators, sources #88, #91 and #97 are Class I 
YSOs; sources #103 and #107 are probably intermediate be- 
tween Class I and Class II; #102 is a Class II; and #92 is inter- 
mediate between Class II and Class III. Magnitudes are missing 
for three sources; we have used the slope of the SED between 
3.6 yum and 24 yum to classify them. Sources #101 and #105 have 
strongly rising SEDs (slopes 1.9 and 1.6 respectively) charac- 
teristic of Class I. Source #1 14 is a Class II source (decreasing 
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1 7^2^50= 



Fig. 13. Condensation 4: The background image is the Spitzer- 
MIPSGAL image at 24 fim; the contours show the 870 /jm emis- 
sion. The contour levels are 0.2 Jy/beam (used to estimate the 
mass of the condensation), 0.4, 0.6, and 0.8 Jy/beam. 




Fig. 14. Condensation 5: The background image is the Spitzer- 
MIPSGAL frame at 24 yum; the contours correspond to the 
870 jjm emission. The contour levels are 0. 1 Jy/beam (used to es- 
timate the mass of the condensation), 0.3, 0.5, and 0.7 Jy/beam; 
the dashed line shows the adopted southern limit of the conden- 
sation. 



slope of -1.15). The nature of #113 is very uncertain, probably 
because it lies too close to the very bright #1 14. 

We have used the Robitaille et al. SED fitting tool to model 
the SEDs of #88 and #105. The parameters of the best fit model 
are given in Table [3] These two sources are rather similar, in 
an evolutionary stage 1, with a massive disk. We tried to model 
#103 and #107 but the results are very uncertain (both stage 1 and 
stage II are obtained); far-IR measurements are clearly needed to 
constrain their SEDs. 

• Condensation 6: This lies on the north-west border of 
RCW 120, apparently far from the main IF. However a north- 
south 24 yum filament, probably tracing the emission of dust in 
the PDR, is seen adjacent to this condensation (Fig. 15 i. Also, 
diffuse Ha emission is observed nearby (Fig.|2]l. Thus this con- 
densation is submitted to the pressure of the ionized gas. Its mass 
is in the range 52-90 Mq; its peak column density corresponds 
to a visual extinction of 18 mag. It appears as an IRDC on both 
the 8.0 fim and 24 fim images. 

Fig. [TS] shows the 870 fj.m emission contours superimposed 
on the 24 //m image. Condensation 6 contains a small cluster of 
at least three YSOs, separated by ~0.17 pc (in projection). YSO 
#106, very bright at 24 //m (possibly saturated), is also detected 
at 70 fim; it lies 0.08 pc away from the dust emission peak. All 
the indicators show that these three YSOs are Class I (#108) or 
flat spectrum sources (#106, #1 10). 

We tried to model the SED of object #106 using the 
Robitaille et al. fitting tool. The evolutionary stage of this source 
is uncertain: the best model, with a central object of mass 2 Mq 
and an accreting envelope, corresponds to stage 1 (Table [3]l, but 
other good models have more massive central objects and no ac- 
creting envelopes, and thus correspond to stage 11 sources. 

The other stars in the field of Fig.[T5]are Class II (#99, #100) 
or intermediate Class Il-Class III (#116). Fig. 15 also shows a 



H = 11.99, K = 11.65) corresponding to a B9V star with a 
visual extinction of the order of 6.1 mag. We suggest that the 
extended 24 /im emission comes from a PDR excited by this star. 
Region C and its exciting star are very similar to regions A and 
B in condensation 4, which we believe to be second-generation 
objects. 




Fig. 15. Condensation 6 and its surrounding. The background 
image is the Spitzer-MIPSGAL frame at 24 //m. The contours 
correspond to the 870 yum emission. The contour levels are 
0.2 Jy/beam (used to estimate the mass of the condensation), 
0.4, 0.6, 0.8, and 1.0 Jy/beam. 



small faint region of 24 fim extended emission which we call 
region C. Its central star has 2MASS magnitudes (J = 12.70, 
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Table 3. Parameters of the "best fit models" for the SEDs of a few sources, obtained using the SED fitting tool of Robitaille et 
al. (120071) 



YSO M„„- 


^star 




^'^envelope 


L comment 


(Mq) 


(K) 


(Mq) 







#38 


3.2 


4350 


2.4 10"' (10^3- 


-10 


') 


1.7 lO-'' 


98 


stage I 


#61 


0.3 


3370 


4.3 10-2 






7.3 10-5 


20 


stage I 


#51 


1.3 


4355 


6.3 10-- 






1.8 10-" (10-*~10-*) 


17 


stage I-Il 


#59 


6.6 


19700 


1.0 10-5 (10-"^- 


-10 


') 





1320 


stage II-III 


#67 


5.1 


7554 


1.7 10-' (10-5- 


-10 


') 


2.2 10-' 


280 


stage II 


#88 


4.1 


4400 


4.4 10-2 (10-'- 


-10 


') 


3 lO-'' 


130 


stage I 


#105 


1.7 


4200 


4.5 10-2 (iQ-3_ 


-10 


2) 


6 


38 


stage I 


#106 


2.0 


4160 


2.210-' 






1.2 10-5 


105 


uncertain 



5.5. The structure of the PDR and the long-distance 
influence of the Hii region 



At first glance RCW 120 appears to be a perfect spherical bub- 
ble. However, as stressed by ZA07, the ionized gas is trying 
to flow out of the bubble, and consequently the Hii region is 
elongated towards the north, where the surrounding interstellar 
medium is of lower density. This morphology is delineated by 
the shape of the ionization front and by that of the low intensity 
870 fim emission zones, as shown in Fig. 
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a composite image 
of the 870 yum emission in turquoise and of an unsharp-masked 
image at 24 yum in red. The unsharp masking emphasizes the 
small-scale structures of the 24 //m emission, both in the ionized 
gas and in the PDR. The black arrows at the top of Fig. 16 point 
to low intensity 870 fj.m emission zones which close the northern 
extension of the ionized bubble; these zones are associated with 
24 jum filaments in the adjacent PDR. The upper aiTow shows 
that a low intensity 870 //m emission region joins condensation 
5 and condensation 9. Here again we are observing the layer of 
neutral material collected around the ionized region. 



Fig. 17 shows the extent of the PDR. This is a colour com- 



posite image of the Ha emission, in red, and of the unsharp- 
masked image at 24 yum in turquoise. We see in red the Ha emis- 
sion of the ionized gas inside the bubble, and in pink a zone of 
diffuse low intensity Ha emission from the PDR surrounding the 
bubble. This ionized PDR is extended. The unsharp-masked im- 
age at 24 fim shows the structures of the dust emission from the 
PDR, covering about the same zone as the diffuse Ha emission. 
The ionized PDR extends over an area with about twice the ra- 
dius of the central Hii region. As stressed by ZA07, the IF is 
porous, and the ionizing radiation escapes to the outside through 
small holes. This ionized PDR is probably of low density, but 
its existence demonstrates the long-distance influence of the H ii 
region and its central exciting star. 

The 870 /urn emission map shows radial features inside con- 
densation 5 and also in the southern parts of RCW 120; there two 
bubbles are observed between condensations 1 and 2 (see the 
black arrows at the bottom of Fig. [T6] l. Fig. |4]shows, at 8.0 fim 
and with a high resolution, the radial structures inside conden- 
sation 5; only the pressure of the ionized gas seems able to form 
filaments perpendicular to the IF. 

Thus star formation can probably be triggered far from the 
ionization front - for example in condensation 6, which is in 
contact with the extended PDR. 




Fig. 16. Colour composite image, with the 870 /im emission in 
turquoise, and in red an unsharp-masked image at 24 fj.m which 
emphasizes the small-scale low intensity structures of the dust 
emission from both the ionized region and the surrounding PDR. 
The black aiTows at the top of the figure point to the collected 
neutral material closing the northern extension of the ionized 
bubble. The arrows at the bottom of the figure point to the bub- 
bles formed by the ionized gas flowing away from the central 
H 11 region. 



5.6. An unanswered question 

A question raised by ZA07 remains unanswered: what is the 
three-dimensional geometry of the RCW 120 bubble? Is there 
a preferential plane for the accumulation of neutral material dur- 
ing the expansion of RCW 120? We have identified only two 
Class I sources in the direction of the ionized gas (#61 and #65), 
whereas more should be observed if star formation is at work 
in a spherical collected shell. According to the simulations of 



14 



L. Deharveng et al.: Star formation around RCW 120, the perfect bubble 




Fig. 17. Colour composite image, with the Ha emission of the 
ionized gas in red, and the unsharp-masked image at 24 fim in 
turquoise. 



Krumholtz et al. (12007b . the presence of a magnetic field in- 
hibits motion of the gas across magnetic field lines, leading to a 
concentration of collected material in a preferential plane. This 
point, important for the study of star formation triggered by H n 
regions, will possibly be addressed via submillimetre polarimet- 
ric observations. 



6. Conclusions 

In this paper we present new observations, obtained with 
the APEX-LABOCA camera, of the thermal dust emission at 
870 fim. The higher angular resolution and signal-to-noise ratio 
of these observations, compared to the previous 1 .2-mm obser- 
vations, allow the detection of low brightness structures and per- 
mit more accurate mass estimates. The observations confimi the 
presence of a shell of dense neutral material collected around 
the ionized region during its expansion. This shell is now ob- 
served all around RCW 120, even at the border of its northern 
extension. The shell's mass is greater than 2000 Mq; we con- 
firm that it is fragmented, with massive fragments of a few tens 
or hundreds of solar masses. We also confirm that a very massive 
(~ 250 Mq) and compact (0.095 pc x 0.05 pc) core lies adjacent 
to the ionization front in the most massive fragment. 

We present and discuss mid-lR observations at 24 jim and 
70 lum, obtained by Spitzer-MIPSGAL. We have measured 138 
sources at 24 fim and six at 70 jum. These observations provide 
new indicators to better distinguish between stars and YSOs, and 
to better distinguish between the evolutionary stages (Class 1 or 
Class II) of the sources. Robitaille et al. (2006) have shown the 
uncertainties tainting these classifications. RCW 120 is a good 



illustration of the difficulty of this study. The main uncertainty 
is due to the objects' external extinction, which is not known; 
highly reddened stars can be mistaken for Class II sources, and 
highly reddened Class 11 sources for Class I sources. We have 
used the 870 jum emission to determine the maximum external 
extinction which may affect a given 24 jum source. 

We have detected 21 Class I sources, 18 sources possibly in- 
termediate between Class I and Class 11, 29 Class II sources, and 
nine more uncertain Class Il-Class III sources. A large fraction 
of the Class I or intermediate Class I-Class II sources are ob- 
served in the direction of the shell of collected neutral material 
surrounding the ionized region. This distribution demonstrates 
that their formation has been triggered by the expanding H ii re- 
gion. 

As stressed by ZA07, we believe that several different mech- 
anisms are at work simultaneously to form stars in the collected 
layer. The mechanisms that we tentatively identify are small- 
scale gravitational instabilities in condensation 1, dynamical in- 
stabilities of the IF in condensation 2, and interaction of the IF 
with a pre-existing molecular condensation in condensation 4. 
The collect & collapse process, which implies large-scale grav- 
itational instabilities of the collected layer along its surface, is 
possibly at work in the dense core of condensation 1, which har- 
bours a source detected only at 70 //m. High-resolution kine- 
matic obsei-vations are necessary to ensure that a Class source 
is present there, building its mass by accreting material in the 
dense core; enough mass is present there (~ 250 Mq) to allow 
the formation of a massive star. 

We have also detected, at 24 jjm, in the direction of con- 
densation 1, a very unusual and remarkable feature: a chain of 
about eleven Class I or flat spectrum sources, parallel to the IF. 
These YSOs have approximately solar masses, and are regularly 
spaced by some 0.1 pc. To our knowledge it is one of the best 
illustrations ever found, in the domain of star formation, of the 
Jeans gravitational instability at work. 

Finally, we have shown that the H n region, via the UV ra- 
diation of its central exciting star, has a long distance influence 
upon its suiTounding medium. If R is the radius of the H ii region, 
the low density ionized FDR, and the associated dust emitting at 
24 fim, extends a further R outward from the IF. This influence 
is also reflected by the shape of the collected layer, which shows 
radial features and bubbles. We conclude that the H ii region can 
probably trigger star formation far from the IF. This remains to 
be confirmed by other examples of H ii regions. 
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